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Abstract 

The branching fractions and the decay amplitudes of B — > J/ipK* decays are 
measured in a 29.4 fb _1 data sample collected with the Belle detector at the 
KEKB electron-positron collider. The decay amplitudes of helicity states of the 
J I ipK* system are determined from the full angular distribution of the final state 
particles in the transversity basis. The branching fractions are measured to be 
(1.29±0.05±0.13) x 1(T 3 for neutral mesons and (1.28±0.07±0.14) x 1(T 3 for charged 
mesons. The measured longitudinal and transverse (perpendicular to the transver- 
sity plane) amplitudes are |A | 2 = 0.62 ± 0.02 ± 0.03 and |A ± | 2 = 0.19 ± 0.02 ± 0.03, 
respectively. The value of |Aj_| 2 shows that the CP even component dominates in 
the B° -► J/^K*°(K s tt°) decay. 

PACS: 13.25.Hq, 14.40.Nd 
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1 Introduction 



The decay B° — > J/ipK* (K* — > Ksir ) has the same quark level diagram as 
B° — > J/ipKs and may, therefore, be used for a determination of the CP vio- 
lation parameter sin 20! . However, since the J/vpK* system has three possible 
helicity states with different superpositions of CP eigenstates, the dilution of 
the CP asymmetry depends on the degree of polarization of the decay. The 
mix of CP even and odd states must be known before a measurement of sin 20i 
can be attempted. 

The CP mix can be determined through a full angular analysis of the final 
state particles to obtain the decay amplitudes of the system in the transversity 
basis. Once the CP mix is measured, the CP eigenstates of the events may 
be projected out in a statistical way. The full angular analysis also serves as 
a test of the factorization hypothesis by measuring the imaginary part of the 
decay amplitudes. 

In this paper, we report the measurements of branching fractions for B° — > 
J/tpK*° and B + — > J/ipK* + and the decay amplitudes of the system ob- 
tained by a full three-dimensional angular analysis. The determination of 
sin20i based on the reconstruction of these decays is described elsewhere[l]. 



2 Data sample 

The data sample used in this analysis corresponds to an integrated luminosity 
of 29.4 fb _1 recorded with the Belle detector[2] at the KEKB electron-positron 
collider [3]. Events are required to satisfy the hadronic event selection criteriafl] 
and have R 2 < 0.5, where R 2 is the ratio of the second to zeroth Fox- Wolfram 
parameters [4]. 

The branching fractions and decay amplitudes are measured by reconstructing 
neutral and charged B mesons in B° — > J/ipK* and B + — > J/ipK* + (inclusion 
of charge conjugate modes is implied throughout this paper). The reconstruc- 
tion of the J/ifi is performed using the dilepton decays, J/vp — > e + e~ and 
fi + [i~ . For the e + e~ mode, electrons and positrons are identified by matching 
between the energy measured in the electromagnetic calorimeter (ECL) and 
the momentum measured in the central drift chamber (CDC), the shape of 
the cluster energy deposit in the ECL, the dE/dx measured in the CDC, and 
the light yield in the aerogel Cerenkov counters (ACC). A likelihood is calcu- 
lated from these measurements and required to be consistent with the electron 
hypothesis. To correct for energy lost by final state radiation, the energy of 
any cluster in the ECL within 0.05 radians of the track momentum is added 
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to that of the track. The invariant mass of each pair of identified electrons is 
calculated, and the pair is identified as J/ip — > e + e~ if the mass is in the range 
2.95 GeV/c 2 < M(e + e~) < 3.15 GeV/c 2 . Tracks are identified as muons by 
means of a likelihood number based on (i) a comparison of the number of layers 
with associated hits in the muon detector (KLM) with the number expected 
based on momentum and (ii) the energy of the associated hit in the ECL. An 
oppositely charged pair of identified muons is identified as J/ip — > [i + fi~ if the 
invariant mass is in the range 3.05 GeV/c 2 < M{jj, + ji~) < 3.15 GeV/c 2 . To 
improve the momentum resolution, a kinematic fit which uses the J/ip mass 
as a constraint is performed on J/tp candidates passing the above selections. 

K*° and K* + candidates are reconstructed in the decay modes K*° — > K + tt~ , 
K*+ -> K s ir + , K*+ -> K+ir°, and K*° -> X 5 7r°. Charged kaons are iden- 
tified by requiring the kaon likelihood of a track to be consistent with ex- 
pectations. The kaon likelihood is obtained by combining measurements of 
the time of flight measured by the scintillation counters (TOF), dE/dx by 
the CDC and the hit information in the ACC. Tracks which are not identi- 
fied as kaons and not used as leptons in the J/ip reconstruction are treated 
as charged pion candidates. Ks candidates are reconstructed from pairs of 
oppositely charged tracks which satisfy three conditions: 1) the distance of 
the closest approach of both tracks to the nominal interaction point is larger 
than 0.03 cm, 2) the angle between the K$ momentum vector and the vec- 
tor displacement of the K$ vertex point from the J/ip vertex is less than 
0.15 radians, and 3) the reconstructed decay vertex of the K s is at least 0.1 
cm away from the interaction point. Each pair with invariant mass satisfying 
0.47 GeV/c 2 < M(n+n-) < 0.52 GeV/c 2 is identified as a K s . To improve 
the mass resolution the momenta are re-fitted, constraining both tracks to 
originate at the reconstructed vertex. Candidate n° mesons are reconstructed 
from clusters in the ECL that are unmatched to charged tracks and have 
energy greater than 40 MeV. A photon pair with an invariant mass in the 
range 0.125 GeV/c 2 < M(jj) < 0.145 GeV/c 2 is identified as a n°. A mass- 
constrained fit is performed to obtain the momentum of the ir°. A K* is 
identified if the absolute difference between the invariant mass of an identified 
Kit pair and the nominal K*(892) mass is less than 75 MeV/c 2 . 



3 Measurement of branching fractions 

Candidate B° and B + mesons are reconstructed by selecting events with a J ftp 
and a K* and examining two quantities in the center-of-mass of the T(4S), 
the beam-constrained mass (M& c ) and the energy difference between the B 
candidate and the beam energy (AE). The beam-constrained mass, which is 
the invariant mass of a reconstructed J / ip and K* calculated taking the energy 
to be the beam energy, is required to be in the range 5.20 — 5.29 GeV/c 2 . For 
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modes with charged tracks only, |Ai?| is required to be less than 30 MeV. For 
decay modes that include a 7r°, AE is required to satisfy —50 MeV < AE < 
30 MeV. To eliminate slow tt° backgrounds, the angle of the kaon with respect 
to the K* direction in the K* rest frame, 9 K *, is required to satisfy cos9 K * < 
0.8. This is equivalent to a requirement that the 7r° momentum be greater 
than 175 MeV/c. When an event contains more than one candidate passing 
the above requirements, the combination for which AE is closest to zero is 
selected. Defining the signal region as 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 , 
we find that the numbers of events passing the selection criteria and lying in 
the signal region are 968 for J/ijK* (K+rr~), 241 for J/iPK* + (K S 7T + ), 220 for 
J/^K* + (K + n°) and 42 for J/^K*°(K S 7r°). 

For each reconstructed mode, the distribution in M bc is binned and fitted to 
a parameterized function of the form 

N(M bc ) = f sig (M bc ) + ]T f l cf (M bc ) + f nr (M bc ) + f combl (M bc ), (1) 



where f sig , f^, f nr , and f CO mU are the shapes for the signal, background from 
other B — > J/ipK* modes (which we will refer to as cross- feed), background 
from non-resonant B — > J/ipKir, and combinatorial background, respectively. 
f S ig is a Gaussian with the peak position fixed at the B mass and a width 
that is determined from Monte Carlo. Only the amplitude is varied in fit- 
ting. f l c f is described by a separate Gaussian for each of the three J/ipK* 
decay modes other than the signal mode. Monte Carlo simulations are used 
to determine the mean and width for each mode, as well as the efficiency 
for it to pass the signal mode criteria. The amplitudes of are scaled such 
that the ratios of cross-feed to signal efficiency are fixed assuming equal pro- 
duction rates and decay branching fractions for the neutral and charged 5's. 
The rate of cross- feed from B — > J/vpK modes is also estimated using Monte 
Carlo and found to be negligible. The total amount of the cross-feed is esti- 
mated to be 1.8-6.2% (depending on the mode) of the events in the signal 
region. The shape of f nr is also assumed to be a Gaussian. While the mean 
and width are expected to be similar to those of f sig , they are determined 
separately to allow for minor kinematic differences, using events in the K* 
mass sideband (1.1 GeV/c 2 < M(Kir) < 1.3 GeV/c 2 ). The magnitude of the 
non-resonant contribution is derived from fitting the Kir mass distribution for 
events selected by the criteria used for the signal reconstruction without requir- 
ing the Kn mass to be inside the K* mass window. The Kir mass distribution 
(Fig. 1) is fitted with P-wave and D-wave Breit-Wigner functions describing 
the K*(892) and .£O!;(1430) mass peaks, respectively, and a threshold function 
constrained at the mass threshold to describe the non-resonant production. 
The peak positions and widths of the Breit-Wigner functions are fixed to the 
nominal PDG values[5]. The amplitude of f nr is determined from the con- 
tribution of the fitted background functions in the K* mass that takes into 
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account the contamination of the combinatorial background in the functions. 
The estimated contribution of the non-resonant decay in the signal region is 
5.2%. f com bi is the line shape for the combinatorial background and is modeled 
by the "ARGUS function" [6] , where the parameters are free in the fit. Fig. 2 
shows the results of the fits to four modes. The fits indicate the contamination 
of the combinatorial background is less than 1% in the signal region. 

The branching fraction is calculated from the yield in the signal region ob- 
tained from the fitted f sig . The detection efficiency for each mode is esti- 
mated by applying the selection criteria to Monte Carlo event samples. Events 
are generated unpolarized in the J/ipK* system using the QQ event gen- 
erator [7] and are passed through a detector simulation program based on 
GEANT3[8]. The estimated efficiencies are 28.4%, 20.7%, 11.8%, and 8.1%, 
for J/ipK*°(K+ir-), J/tpK* + (K s ir + ), J/4>K* + (K + n°) and J '/ > K*° (K s n ) , 
respectively. The branching fractions of secondary decays are fixed to the 
PDG values [5] and the numbers of neutral and charged BB pairs from T(4S f ) 
decays are assumed to be equal. 

The systematic uncertainties on the branching fractions are dominated by the 
following sources: 1) tracking efficiency (2.0% per track), 2) n° reconstruc- 
tion (3.0%), 3) lepton identification efficiency (5.0%), 4) kaon identification 
efficiency (< 1.0%), 5) background estimation (2.3-4.8% depending on the 
mode), 6) number of BB events (1.0%), 7) polarization (1.7%), and 8) branch- 
ing fraction for secondary decays (1.7%). The estimate of the non-resonant 
B — > J I if)K7r background has a relatively large uncertainty because the fit to 
the Kn mass distribution is poor in the region of Kit mass from 1.1 GeV/c 
to 1.3 GeV/c, where some excess is observed. This excess has also been seen 
elsewhere [9, 10], and there has been some speculation as to its origin[9]. As its 
source is not established, we estimate the uncertainty in its contribution to 
the signal by including in the fit to the Ktt mass distribution an additional 
term that is constrained to zero at the mass threshold and extends linearly to 
1.35 GeV/c 2 . The change to the signal yield due to this modification is taken 
to be the uncertainty. 

The resulting branching fractions are summarized in Table 1 with the esti- 
mated systematic errors. 



4 Measurement of decay amplitudes 

The decay amplitudes of B — > J / ipK* decays are measured in the transversity 
basisfll]. The direction of motion of the J/ip in the rest frame of the B can- 
didate is defined to be the x-axis. The y-axis is defined to be in the direction 
of the projection of the K momentum into the plane perpendicular to the 
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x-axis in the B rest frame. The x-y plane contains the momenta of the J/ip, 
the K, and the ir. The z-axis is defined to be perpendicular to the x-y plane 
according to the right-hand rule. The angle between the / + direction and the 
z-axis in the J /if) rest frame is defined as 9 tr . The angle between the x-axis 
and the projection of the l + momentum onto the x-y plane is defined as (p tr in 
the same frame. The angle 9k* is defined as described in the previous section. 

The distribution of these three angles for B — > J/ipK* decays is described in 
terms of three amplitudes [12]: 



1 dT 9 

rl a— a a a a = x [2 cos 2 (1 - sin 2 ^ r cos 2 tr ) | A | 2 

T a cos 9 tr d cos 9k* d<pt r 327r 

+ sin 2 9 K * (1 - sin 2 9 tr sin 2 tr ) | A\\ | 2 

+ sin 2 9 K * sin 2 ^ r |Aj_| 2 

+i] sin 2 6k* sin 29 tr sin (j) tr Im(A^A±) 

— \= sin 29 K * sin 2 9 tr sin 2(j) tr Re(AlA\\) 
v2 

+rj—j= sm29 K * sm29 tr cos (j) tr I m(AlA±)], (2) 
v2 



where Aq, A\\ and A± are the complex amplitudes of the three helicity states 
in the trasversity basis, and rj = +1 (—1) for B° and B + (B° and B~). 
\Aq\ 2 denotes the longitudinal polarization of J/ip while |Aj_| 2 (|A|| | 2 ) is the 
transverse polarization component along the z-axis (y-axis). Also, |A | 2 + |v4||| 2 
is the amplitude corresponding to the CP-even state, while |Aj_| 2 is the CP-odd 
component in B° — > J/v/jK*°(K*° — > K s tt°). These amplitudes are normalized 
so that: 

|A | 2 + |A||| 2 + |A ± | 2 = 1. (3) 



The amplitudes are determined by fitting this function to the measured three- 
dimensional distribution in 9 tr , <p tr and 9k*, taking into account the detection 
efficiency and background. The resolution of the angular measurements is 
estimated by Monte Carlo and found to be typically less than 0.02 radians. 
The value of 77 is determined from the charge of pions or kaons used in the K* 
reconstruction. We do not include the B° — > J/^K*°(K s tc ) mode in the fit 
since in this case r] is not well-defined. 

The fit is performed using an unbinned maximum likelihood method. The 
probability density function (PDF) is defined using the theoretical distribution 
in (2) and can be expressed as 
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PDF(x, y, z, M bc ) = N x [f sig {M bc ) x e(x, y, z) x - ——— (a;, y, z) 

T dxdydz 

+ J2fc f (M bc )xADF cf (x,y,z) 

i 

+ fnr(M bc ) x ADF nr (x,y,z) 

+ fcombi(M bc ) x ADF combi (x,y,z)], (4) 

where x = cos6 tr , y = <ptr, and z = cos6k*, N is the normalization factor of 
the PDF, e(x, y, z) is the detection efficiency as a function of the three angles, 
and ADF c f, ADF nr , ADF combi are the angular shapes for the cross- feed, non- 
resonant and combinatorial backgrounds, respectively. 

Here, f sig {M bc ), f l cf (M bc ), f nr (M bc ), and f com bi(M bc ) are the fractions of signal 
events, cross-feed contamination, non-resonant contamination and combina- 
torial background as a function of beam constrained mass, respectively. These 
fractions are obtained from the measurements of the branching fractions de- 
scribed in the previous section. 

The detection efficiency function e(x, y, z) is obtained from a large Monte Carlo 
sample of 1 million events for each mode generated without any polarization 
of the B — > J/ipK* system. Events are histogrammed in a 20 x 20 x 20 grid in 
a cos9 tr — (ptr — cosOk* cube. The distribution is fitted to a three-dimensional 
polynomial with correlations taken into account. The efficiency is almost flat 
except in the region cos9 K * ~ 1) where the pion is slow so that the efficiency 
is reduced. 

The angular distribution function for the cross-feed background (ADF c f) is 
determined from the distribution in Monte Carlo events. The function for the 
non-resonant production (ADF nr ) is determined from events in the sideband of 
the Kn mass distribution. The distribution for the combinatorial background 
(ADF combi ) is obtained from events in the sideband of the beam-constrained 
mass. These distributions are parameterized as polynomials where the param- 
eters are determined from the fit. 

In the fit to the angular distributions, the imaginary part of A is defined to 
be zero relative to the imaginary parts of the other amplitudes since the over- 
all phase of the decay amplitudes is arbitrary. By applying the normalization 
condition (3), four parameters, |A)| 2 , |^4-l| 2 > ar 9(An) and arg(A±) are left be 
determined from the fit. The normalization of the PDF (N) is calculated by 
numerical integration over the three dimensional angular cube whenever pa- 
rameter values are changed. A likelihood is defined as the product of the PDF 
for all the events and the values of parameters are determined by maximizing 
the likelihood. We perform separate fits to the three decay modes as well as a 
combined fit by defining a single likelihood. The parameter values determined 
from these fits are summarized in Table 2. The projected angular distributions 
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of the data for all K* modes combined are shown in Fig. 3. The distributions 
are corrected for the effects of detector acceptance and backgrounds. 

Systematics in the fit are studied for the following uncertainties: 1) efficiency 
function (Monte Carlo statistics and effect of polarization), 2) angular dis- 
tribution functions for backgrounds, 3) background fractions, 4) slow pion 
efficiency, and 5) fit algorithm. These contributions to the systematic error 
are summarized in Table 3. The dominant contributions arise from the de- 
pendence of the efficiency and background functions on the polarization, and 
also the uncertainty in the detection efficiency for slow pions. The first un- 
certainty is estimated by comparing the functions obtained for Monte Carlo 
samples generated with and without the polarization in the J/ipK* system. 
The second uncertainty dominates in the cos 9k* distribution where cos 9 k* is 
close to 1. The uncertainty is estimated by comparing with the results with a 
cut at cos^x* < 0.9 for decay modes with n^ 1 and with a cut at cos 9k* < 0.7 
for the mode with ir°. The systematics in the fit algorithm are studied using 
toy Monte Carlo events. The differences between input decay amplitudes for 
Monte Carlo and the values obtained by the fit are taken as the uncertainties. 



5 Conclusion 

Branching fractions for B° — > J/ipK*° and B + — > J/ipK* + are measured 
via reconstruction of K*° -> K + ir~, K*° -> K s tt°, K* + -> K + tt°, and 
K* + — > K s ix + with leptonic decays of J/ip. The resultant branching frac- 
tions are Br(B° -> J/^if* ) = (1.29 ± 0.05 ± 0.13) x 10~ 3 and Br{B+ -> 
J/i/jK* + ) = (1.28 ± 0.07 ± 0.14) x 10~ 3 . These values are consistent with the 
measurements by CLEO[10]and BaBar[13] (Table 4). 

The decay amplitudes for B — > J/ipK* are measured by fitting the angular 
distribution of final state particles in the transversity basis to a theoretical 
distribution with the effects of the detector acceptance and the background 
taken into account. In the fit, the sum of amplitudes is normalized to 1 (|A | 2 + 
|^4||| 2 + |^4_l| 2 = 1) an d arg(A ) is defined to be zero. The fit is performed using 
an unbinned maximum likelihood method. From the fit, the values obtained 
are \A \ 2 = 0.618 ± 0.020 ± 0.027, \A ± \ 2 = 0.191 ± 0.023 ± 0.026, arg(A\\) = 
2.83 ± 0.19 ± 0.08, and arg(A ± ) = -0.09 ± 0.13 ± 0.06. The measured value 
of | ^4 j_ 1 2 shows that the CP even component dominates in the decay B° — > 
J/i/jK* (K S Tc ). The value is used in the determination of sin20x as described 
in [1] . Table 5 shows the comparison with other measurements. The parameter 
arg(A\\) is sensitive to final state interactions (FSI). A shift from or n shows 
the existence of the FSI. The measured value is consistent with results from 
other experiments [9] [14]. However, the significance of its deviation from n is 
not sufficient to conclude the existence of FSI. 
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Table 1 

Measured branching fractions for B° — > J/t/jK*° and 1? + — ► J/ipK* 



Mode 


Branching Fraction (x 10 3 ) 


B° -► J/iPK*°(K + ir-) 
B° J/i/jK*°(K s tt ) 


1.30 ±0.05 ±0.12 
1.20 ±0.20 ±0.15 


B° -► J/^if* 


1.29 ±0.05 ±0.13 


J/ipK* + (K s ir+) 
B + -► J/i>K* + (K + ir°) 


1.23 ±0.11 ±0.15 
1.33 ±0.10 ±0.14 


B+ -► J/ipK*+ 


1.28 ±0.07 ±0.14 



Table 2 

Measured decay amplitudes for £>° — ► J/if)K*° and .B + — > J/^K* + (statistical error 
only). 



Mode 


|4>| 2 




ar(/(A||) 




5° -► J/^K*°(K + tt-) 
B+ -► J/%bK* + (K s ir+) 
B+ -► J/ijK* + (K + TT°) 


0.607 ±0.020 
0.630 ± 0.048 
0.653 ±0.052 


0.195 ±0.023 
0.171 ±0.057 
0.194 ±0.087 


2.87 ±0.19 
2.74 ±0.39 
2.69 ± 1.08 


0.04 ±0.13 
-0.45 ±0.33 
-0.23 ±0.56 


Combined 


0.617 ±0.020 


0.192 ±0.023 


2.83 ±0.19 


-0.09 ±0.13 



Table 3 

Systematic errors in the measurement of decay amplitudes. 



Item 


|4>| 2 


l^l 2 


arg(A\\) 


arg(A±) 


Efficiency 


0.008 


0.003 


0.04 


0.03 


ADF for backgrounds 


0.007 


0.003 


0.01 


0.04 


Background fractions 


0.001 


0.001 


0.01 


0.01 


Slow pion efficiency 


0.025 


0.026 


0.07 


0.02 


Fit algorithm 


0.001 


0.001 


<0.01 


0.01 


Total 


0.027 


0.026 


0.08 


0.06 



Table 4 



Comparison with previous measurements of branching fractions for B° — > J/tpK*' 
and B + — > J/ipK* + . The first (second) quoted error is statistical (systematic). 





5r(£° -► J/i>K*°) (xl0~ 3 ) 


Br(B+ -► J/i>K*+) (xl0~ 3 ) 


CLEO[10] 
BaBar[13] 


1.32 ±0.15 ±0.17 
1.24 ±0.05 ±0.09 


1.41 ±0.20 ±0.24 
1.37 ±0.09 ±0.11 


This measurement 


1.29 ±0.05 ±0.13 


1.28 ±0.07 ±0.14 
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Table 5 

Comparison with previous measurements of decay amplitudes in B — > J/ipK* . The 



first (second) quoted error is statistical (systematic). 



Group 


IA)I 2 


l^l 2 


arg(A\\) 


arg(A±) 


CLEO[10] 
CDF [14] 
BaBar[9] 


0.52 ±0.07 ±0.04 
0.59 ±0.06 ±0.01 
0.60 ±0.03 ±0.02 


0.16 ±0.08 ±0.04 

0.13^ ±0.06 
0.16 ±0.03 ±0.01 


3.00 ±0.37 ±0.04 

2.2 ±0.5 ±0.1 
2.50 ±0.20 ±0.08 


-0.11 ±0.46 ±0.03 

-0.6 ±0.5 ±0.1 
-0.17 ±0.16 ±0.07 


This measurement 


0.62 ± 0.02 ± 0.03 


0.19 ±0.02 ±0.03 


2.83 ±0.19 ±0.08 


-0.09 ±0.13 ±0.06 
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0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 

M(K + n) (GeV/c) 

Fig. 1. The invariant Kir mass distribution for B° — > J/ipK*°(K + Tr~) candidates. 
The solid line shows a fit to two Breit-Wigner functions corresponding to i^"* (892) 
and if|( 1430) with a background function (dashed line). 
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Fig. 2. The distributions of the beam-constrained mass for four reconstructed modes. 
The solid lines show the fits. 



16 






-0.5 0.5 

cos Qtr 



Fig. 3. Projected distributions of cos6 tr , <p tr and cos 9k* for the combined data. 
The distributions are background subtracted and acceptance corrected. The solid 
line shows the theoretical predictions with the obtained decay amplitudes. 



17 



